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Abstract 

This paper proposes a variational principle for the solutions of quantum 
field theories in which the "trial functions" are chosen from the algebra of 
asymptotic fields, and illustrates this variational principle in simple cases. 

1 Introduction 

The most-used methods to find approximate solutions of quantum field theories are 
based on path integrals |l|, H, ||, ^. They have many advantages; however the 
Hilbert space and particle structure of field theory are not evident from this point 
of view. Fock space methods, such as the Tamm-Dancoff approximation]^, ^ and 
the discretized light cone quantization approximation (DLCQ)0 place the Hilbert 
space and particle structure of the theory in the forefront; however the covariance 
of the theory is not evident. The functional Schrodinger picture allows intuitive 
guesses about the form of the solutions of a field theory to be incorporated, but this 
method also fails to be explicitly covariant. 

Another method, less developed than those just mentioned, is the expansion 
in normal-ordered asymptotic fields, the "Haag expansion" or "N-quantum approx- 
imation," which applies directly only to theories without zero-mass particles. In 
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particular, it does not apply directly to gauge theories. Nonetheless, this method 
can be extended to gauge theories; indeed it has been applied to quantum electro- 
dynamics. Using the Haag expansion, one works in the algebra of asymptotic fields, 
and can keep creation and annihilation parts of operators on the same footing, since 
one can choose not to apply elements of the algebra to the vacuum, where terms 
with annihilation operators would annihilate the vacuum and be lost. By contrast, 
in the Tamm-Dancoff approximation, the annihilation parts destroy the vacuum and 
disappear from the calculation. In addition to losing explicit Lorentz invariance, this 
asymmetric treatment of the annihilation and creation parts of the fields destroys 
crossing symmetry. The DLCQ method also treats annihilation and creation oper- 
ators asymmetrically; in addition, it suffers from a lack of explicit covariance. This 
complicates renormalization considerably. 

Just as variational methods have been used in other approaches to quantum 
field theory, this paper proposes a variational principle based on the Haag expan- 
sion. Variational principles in quantum mechanics are powerful ways to go beyond 
perturbation theory. In quantum mechanics, the solutions lie in a Hilbert space 
of functions and the trial functions are chosen from this space. For example, to 
approximate the ground state, one can choose a wave function ip that depends on 
some parameters and determine the parameters by minimizing the ground state 
energy [P]. The exact ground state energy will be less than the approximate energy 
at the minimum. Many attempts have been made to carry over this approach to 
quantum field theories. In field theory, the ground state wave function is replaced 
by a vacuum wave functional that can depend on functions as well as parameters. 
Minimizing the vacuum energy determines the functions and parameters and yields 
approximate information about the solution of the theory |1I0|, pJ] , p^ . This ap- 
proach is usually restricted to wave functionals closely related to Gaussians, because 
the necessary path integrals can be done only in that case. Also, this approach does 
not take advantage of the particle spectrum that we expect to occur in quantum field 
theories. For a relativistic theory, for example, the spectrum should consist of a vac- 
uum of energy-momentum zero, one or more single-particle states (including possible 
bound states) with various masses (here, as in the rest of this paper, the analysis is 
restricted to cases without massless states), and many-particle states whose energies 
and momenta correspond to several massive particles. This particle structure can 



be put into a variational calculation at the outset by choosing to approximate the 
fields, rather than the states, and by using an expansion in asymptotic fields for 
the interacting fields. (Alternatively, as given below, one can use generalized free 
fields in the expansion of the interacting fields and derive the fact that the fields 
in the expansion are ordinary free fields.) Thus, in quantum field theory, assuming 
completeness and irreducibility of the algebra of asymptotic fields, one finds that 
the solutions lie in this algebra and the trial wave functions are replaced by trial 



operators chosen from this algebra. Rudolf Haag |Tj] introduced the idea to use an 
element in the algebra of asymptotic fields to represent the interacting field. It is 
fitting to call such a representation a "Haag expansion." Some applications of Haag 
expansions are given in |jl5|]. While a systematic approximation using terms with 
normal-ordered products having arbitrarily high degree in asymptotic fields leads to 
amplitudes with arbitrarily many momenta and becomes intractable, a variational 
trial operator in the algebra of asymptotic fields can have infinite degree, but still 
can be parametrized in a tractable way. Section 2 formulates the variational princi- 
ple and proves that the minimum picks out the solution of the field theory. Section 
3 illustrate the principle with the simplest free neutral scalar field of mass 

m with the trial operator chosen to be a generalized free field, and a free Dirac 
field where the trial operator is a free Dirac field. To compare with the variational 
method using a vacuum functional. Sec. 4 studies the (j)^ model in 1-1-1 and compares 



the results in the lowest approximation with the results of |jTO[- To show how the 
principle works in cases where the trial element has infinite degree. Sec. 5 applies the 
principle to the gradient coupling model, which has a "nucleon" Dirac field coupled 
to a scalar "meson" field. In that case, one can choose the trial operator for the 
nucleon field to be the product of an arbitrary function of a free scalar field and a 
free nucleon field all at a single spacetime point and the trial operator for the meson 
field to be a free scalar field. Section concludes with the outlook for future work. 

2 The variational principle 

Using reasonable assumptions (in particular that no massless fields are present) and 



a physicist's level of mathematical rigor, one can show||16|| that the asymptotic fields 
obey the same Poincare transformation law as the interacting fields. For example, 
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for a scalar field, 

U{a, A)(j)'''{x)U{a, A)^ = 0"^(Ax + a), (1) 

A G 50(1, 3), a G R^. For the general case, 

U{a, A)ij^'{x)U{a, A)^ = Vij{A-^)ij]'{Ax + a), (2) 

where A in SL{2, C) is replaced by A in 50(3, 1) for integer spin fields. One can also 
show that the asymptotic fields obey free field commutation relations, for example, 
for a scalar field, 

[r{x),rm- = ^H^-y;^'). (3) 

where 

iA(x, m^) = (27r)-^ J d'^ke{k^)6{k^ - m^)exp{-ik ■ x) (4) 

is the Pauli- Jordan commutator function. These two conditions, together with the 
requirement that the vacuum have zero energy, imply that the generators of the 
Poincare group are the free bilinear functionals of the asymptotic fields; in particular 
that (apart from a constant which is the vacuum matrix element of the Hamiltonian) 
the Hamiltonian is diagonalized by the asymptotic fields, 

P'' = H = C0nst.+J2 Hfree W] ■ (5) 
i 

To see this, assume H has an arbitrary expansion in asymptotic fields, 

n n 

H = F(°)+^F!')_,.Af(0) + ^Fj'^'^(0)t 
j=i j=i 

+ E[/ M(k)^r(k)Ar(-k)ea;p(-2^^,x°) ,,(k)Ar(k)tAr(k) 

s,t 

oo „ 73n L 

+fS ,,(k)Ar(k)tAr(-k)texp(2^.;,x°)] + ^ / ^F^''\h) 

X ■.llA'pik.f ■.6{Y,±ki)exp{i±Uk,t±ki-^,), Uk = Vk^ + m^ (6) 

j 

A^'^{ki)^^ stands for either the creation or the annihilation operator, normalized 
relativistically. From Eq.(|^) and the Fourier transform 

0*"(x) = (27r)-=^/2 J d^k4>'''{k)S^{k)exp{-ik ■ x) (7) 
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(j)^'''\k)5m{k) = e{k')—^5{k^ - uj,) + g(-fc°) ^ ^ 5{k' + ^fc), (8) 
the commutation relations of the creation and annihilation operators are 

[A'"(k), t(l)]_ = 2Vk2 + m25(k - 1) (9) 

The nonrelativistically normalized annihilation and creation operators are a*-*"''(kj)*^) = 
(2£^fc)~^/2y4*"(kj)(\ = > 0. Put the form for H into the infinitesimal form of 
Eq.(i, 

i[H,r\x)]=d^r\x). (10) 

The commutation relation, Eq.(^or H), say that a normal-ordered term in H with n 
factors of cj)""'^ will contribute to a term with n — 1 factors of 0"'* on the right-hand- 
side of Eq.(|T^). Since the right-hand-side of Eq.(|TUp is linear in 0"*, the only terms 
allowed in H are those with n = or n = 2. The n = term is the vacuum energy, 
and because 0"* is a free field, the n = 2 term is the free field Hamiltonian. 

The Haag expansion of the interacting fields, stated for simplicity for a single 
neutral scalar field with only scalar bound states, is 



oo -| „ n 

EryE / Y.d'''xj'^-\{x-x.}) : nC(^^) ^ (11) 



{X} = 

n=0 '"• u " i i=l 



where the 0*" include in fields for stable bound states, if there are any. The term 
= t) is a constant that occurs for a scalar field when symmetry is broken. The 
term with /^^^ = 1 is just the in field with coefficient one to fix field strength 
renormalization. In momentum space, the expansion is 



~m = i^nf^'vSik) + 4>'''{k)Srn{k) + E - E / /^"H^l' ^n) 

n n 

X : l[4>Tik^)Sm,ik.^ ■■ 5ik -J2ki)l[d%, (12) 



/(")(xi,...,x„) = (271)-^/'-'^/^ J f('^\k^, ...,kn)expi-tj:ki ■ Xi)Ud%. Hermiticity 
of the field implies f^^^ki, kn) = /(")*(-A;i, -kn). Similar expansions hold for 
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out fields. When this expansion is inserted into the Hamiltonian, the result is an 
infinite series of the form already given in Eq.(^). 

For a given Haag expansion, Eg. (p^ , parametrized by v and the /"-"^'s, the 
F(")'s of Eq.(D are functionals of v and the /"^"^'s. From the discussion above, for 
the exact solution F^"^ = for n = 1 and for all n > 2, F^"^^ = for the and 

(2) 

++ cases, and F^_ is minimum. In other words, as just discussed, the Hamiltonian 
is the sum of free-field Hamiltonians for each in field, together with a constant 
term which is the vacuum matrix element of H. The condition that the Haag 
expansion diagonalizes H leads to an infinite set of nonlinear integral equations in v 
and the f^^^^s. A solution to this set of equations is equivalent to the solution of the 
field theory in all sectors. In practice, it will be difficult to find an exact solution. 
Variational methods can yield approximate solutions. Many different conditions 
can be imposed on the F*^"^'s to find an approximation solution from a variational 
principle. One such condition is to minimize the sum Q (with possible weighting 
factors A„) of the integral of the absolute squares of the coefficients of the operator 
terms, 

Q = Ao|F(°f + A_|f1')|2 + A+|F«|2 

+^—J d'k\FL^l{k,-k)\^ + X+. J d^k\Fi^l{k,k)\^ 

+A++ / d'k\Fi'lik,-k)\^ + J2>^n [ lld'h\F^''\k,,...,k^)\'' (13) 

(here, to simplify the notation, I dropped the subscripts that label the possibly 
different asymptotic fields). The A's are arbitrary positive or vanishing numbers that 
can be chosen to control the weight attached to each term. If the Haag expansion 
of the interacting field has finite degree, the terms in the Hamiltonian of highest 
degrees in in-fields cannot possibly vanish, so the A's for such terms should be 
chosen to vanish. For infinite degree Haag expansions which remove the restriction 
to weak coupling one can keep all the A's positive. For the exact solution all off- 
diagonal terms vanish and the diagonal terms are minimum, so this principle gives 
the exact solution of the theory in all sectors at the absolute minimum if all the A's 
are positive. If all A's except the one multiplying the a^a term for a bound state are 
chosen to vanish, this variational principle reduces to the quantum mechanical one 
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for the bound state. Another possible condition is to minimize the integral of the 
absolute square of the coefficients, 

Jl[d%\F^-\h,...,K)\', (14) 

for some set of values of n and for each combination of creation and annihilation 
parts of the operators. One can also choose to require that some set of sums of the 
positive terms in Eq. ([I3D vanish. Each of these possibilities gives a set of equations 
for the solution of the variational principle. The specific conditions one should 
impose in a given problem should be chosen by experience. 



3 Application to the free field 

The simplest case on which to test this variational principle is the free field. The 
Hamiltonian for a free neutral scalar field is 

H=^J d?x{^'' + (V0)2 + m202) (15) 

As a trial operator, choose a generalized free field, whose two-point function 

has an unknown positive measure p(k^). To avoid the trivial case where p = 0, 
require /g°° p(K^)(iK^ to have a fixed positive value; the exact value doesn't matter. 
Let the field (pgff be represented as the sum of a term, (pdisi^; /i^), with a discrete 
weight at mass p^ and a term with a continuous weight cr, (pconi^', cr), 

(Pgffi^) = <Pdis{x] p^) + (pconix] cr), (16) 
1 r d^h 

<PdUx; f^') = J ^^^[A{k)e-^'- + A\k)e^''-% = E,ik), (17) 

0con(x;a) = (^/^«^'^£|^[^(k;'^')e-^^-^ + S^(k;«:^)e^^-1, 

fco = E^{k) = Vk2 + ^2. (18) 

The relativistically normalized commutation relations are 

[A{k),A\i)]. = 2E,{k)S{k-e), (19) 
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[B(k; K^), B\i; X^)]_ = 2E^(k)S(k - i)6{K^ - \^)a{K^), (20) 

other commutators vanish. To find the minimum of the energy of a given particle, 
one should minimize the coefiicient of the a^a term. The relations between the 
relativistically (capital letters) and nonrelativistically (lower-case letters) normalized 
operators are 



A{k) = J2E,A,\^)a{k), (21) 



B{k;K')^^/2Ejk)b{k-K'). (22) 

When the trial operator (pgff is inserted into the free scalar Hamiltonian, the result 
is 

+2(£;^(k)2 + + m^) : at(k)a(k) :} (24) 



H, 



= / [(-E^(k)E,(k) + k^ + m^) : a(k)6(-k; A^) : 

Xg-i(i?^(k)+E,(k))a:0 ^ [E^{\i)E),{\i) + k^ + m^) 

X : 6t(k; A')a(k) : g-^^^''^-^^^)^" + h.c], (25) 



2-^ 2 J 2^E,{k)Ex{k) 

x{[(-£;,(k)^,(k) + k^ + m') : b{k- K')b{-k- X') : 

^g-i(E„(k)+BA(k))xO ^ ] ^ 2(^,(k)EA(k) + k^ + m^) 

X : 6t(k; A2)6(k; k^) : e-(^^«-^A (26) 

For the terms of the form 1 (vacuum energy) or a^a (particle energy), we must 
minimize the energy (2k^ + //^ + m^) / with respect to /x^, or, alternatively, 
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minimize the same energy written as [i?^(k)^ + + m^]/£'^(k) with respect to 
E^(k); for terms of the form aa or a^a^ we minimize the absolute value squared of 
the coefficient (since we want to bring the Hamiltonian to diagonal form). For terms 
of the form b'^b, we must minimize the energy, 

with respect to and A^. For terms of the form or a^b, similar expressions 
have their minima for a concentrated at = /i^ = m^. For terms of the forms aa, 
a^a\ ab, a^b\ bb, and 6^6^, the squares of the coefficients vanish for fi"^ = and a 
concentrated at = m^. For the contributions from (peon, the minimum occurs for 
cr at = m^. The net result is that the minimum of the operator H occurs at 

(Pgffi^) = <Pdis{x] m^),(j)con = 0; (28) 

i.e., for a free field of the mass in the Lagrangian. This is no surprise. In general, 
for theories without massless fields or particles, we would assume free field form for 
the asymptotic fields without doing a calculation. (Note that these results require 
minimizing an expression that has dimensions of energy.) 

If one assumes a free field in the corresponding calculation of the Hamiltonian 
for the Dirac field the b^b and Sd terms are diagonal in helicity; the bd and b'^d^ 
terms are not. The minimum of the absolute valued squared coefficients also yields 
the result /z^ = m^. 



4 Application to the 0^ model in 1+1 

A less trivial, but still elementary, example is the 0^ theory in one space, one 
time dimension. This model was studied recently using a variational method by 
G. Tiktopoulos ||T0[. (The literature on variational calculations can be traced from 
Pm , |TT], |T^, |T3|.) The example just below illustrates the variational principle in the 
algebra of asymptotic fields for this theory in the lowest non-trivial approximation, 
where the calculations can be done easily by hand. Higher approximations require 



symbolic manipulation programs that have been developed using Reduce 3.5. We 
will report on the results of higher approximations separately. The Lagrangian is 

J!:^lid,<f>-d^<f>-ml<p')-^cf>'. (29) 

The Hamiltonian density is 

n^l{^' + <P" + ml<l>') + ^4>'- (30) 

For a theory that has no massless particles or fields, the (on-shell) asymptotic fields 
will be an irreducible set of operators in which to expand the interacting field. 
The demonstration given above that the free asymptotic fields diagonahze the free 
Hamiltonian supports this expectation. The lowest variational ansatz is then 

(f){x) ^v + (f)^'''\x;m'^), (31) 

where the constant f is a symmetry-breaking vacuum matrix element of that 
may occur due to radiativcly-induced symmetry breaking and 0*^*"-' is a free field of 
unknown mass m. Substituting this ansatz and re-normal ordering, the Hamiltonian 
becomes 

n = [\{<<Pl> + <<i^'o'>)+ml{v'+<<l>l>) 
+ ^{v' + 6v'<cl>l>+3<cl>l >')] 
+[{mlv + ^v^) + < 4 >] : (1)0 : 

+ \ : 0^ +^ : 0^ :, (32) 

where 0o stands for 0(*")(O; w?). The vacuum matrix element is just the first square 
bracket above. The coefficient of the term in (1/2) : 0q : should be the square of the 
physical mass m^, i.e. 

ml^ml + ^{v^+<ct>l>). (33) 
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This is the gap equation in this approximation. Following Tiktopoulos, one can 
verify this below by minimizing the vacuum energy density which is a functional 
of the energy, E{k), of the trial in-field quanta of momenta k with respect to this 
energy. Regulate the vacuum and other divergent expressions with a momentum- 
space cut-off; the two-point function is then 

Wo(.)^o(y)>^^£^e-W--'. (34) 

Then 

The vacuum energy density is 

X r. dk X . dk ... / n 

The minimum is given by 

6{H) _ 1 k'^ + ml X X dk' 1 

SE{k) ~ 8^^ E{ky ^~ 167rE{ky~¥^ J-AE{k')E{ky 

= 0. (37) 

After multiplying by STrE^k)"^, the result is the gap equation, so m? = ml as given 
by (H). 

Since the gap equation relates the physical mass ml, which should be finite, to 
the bare mass and the divergent integrals cut-off at A, one can follow Tiktopoulos 
in introducing another finite mass, M, via 

M' = ml + — f"^ (38) 

This form of renormalization can replace the usual one, since here only the mass is 
renormalized. As Tiktopoulos showed, the gap equation then relates finite quanti- 
ties, 

ml = M' + + ^In^. (39) 
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The minimum of the vacuum energy density is 



(n) ^ — dk k / f + —ml 

A 4 1 2,,, 2 A , M\ 1 . A , M\2 

+ 7T7^ + ^ + ^^^^ + TTT ^ + ^^^^ > 40 
24 2 ^ Stt ml ^ 2A Stt m? ^ ' ^ ^ 



again in agreement with Tiktopoulos. Thus the present variational calculation leads 
to the same results as does that of Tiktopoulos in the lowest approximation. This 
agreement will not persist in higher approximations. Tiktopoulos adds Gaussians 
to his vacuum functional in order to improve his approximation. The Haag expan- 
sion suggests a different form of the higher-order terms, with higher-degree normal- 
ordered products of asymptotic fields, in the expansion of the interacting field. These 
are very different approximations. If the in (or out) fields are irreducible, then the 
Haag expansion will approximate the exact solution in theories without massless 
fields. In higher approximations, the present variational method will place emphasis 
on minimizing the absolute square of the coefficients of the non-quadratic terms 
in Ti and on minimizing the coefficients of the a^a terms which correspond to the 
energies of the asymptotic quanta. 



5 Application to the derivative coupling model in 



To illustrate a trial function that is of infinite degree, consider the derivative coupling 
model, with 



1+3 




(41) 



where (f) and ip are renormalized fields. The Hamiltonian is 




Assume = (pQ, ip =: /(</>o) : 'ipo- The Hamiltonian becomes 



H ^ I d'x[Z,i;o : fM : (^ : f'{Wd^cj>o : +57^5^00 : /(0o) : +i : f{<f>o) : ^'d^ + 
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The first two terms in the bracket will cancel if 



tf + 9f = 0, (44) 

where the fact that {d^cfxp) = allows removing (9-^00 from the normal-ordered 
product in the first term. Field strength renormalization of ip requires /(O) = 1. 
The solution of this is f{x) = exp{igx). Evaluation of : /{(poY :: /(0o) : using the 
Baker-Hausdorff-Campbell theorem gives exp{g'^ {(j){x)(j){y))) in the limit x —>■ y for 
this product . This determines 

^ = lim^^yexp{g^{^{x)(f){y))). (45) 

Thus the solution that reduces the Hamiltonian to free field form is 

ip{x) =: exp{ig(f)o{x)) : ipoix), 4>{x) = (po{x). (46) 

(Antisymmetrization of the Fermi fields ip and ip"^ above complicates the calculation, 
but does not change the outcome.) 



6 Outlook for future work 

Much must be done to make this idea into a useful tool for field theory calculations. 
The method should be applied to examples with a finite number of terms in the 
Haag expansion that serves as the trial element. Calculations of vertex functions in 
an all-scalar relativistic model with the interaction are ongoing. This method 
will be used to calculate hydrogen bound states with the soft photon cloud taken 
into account for the charged particles. The full power of the method will become 
apparent only when infinite degree Haag expansions, parametrized in a tractable 
way, are used. The derivative coupling model discussed above is a trivial example 
of such a parametrization. As mentioned above, the soft photon cloud around a 
charged particle will be represented using asymptotic fields in a later paper. Realistic 
models will be much more difficult to treat. We must also confront the problem of 
the coupling of high-energy and low-energy modes pointed out by Feynman||17||. The 



step to nonabelian gauge theories in which the confined fields do not have asymptotic 
fields will be the most difficult step. Hopefully the case of electrodynamics, where 
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the charged fields acquire a cloud of soft photons, which we believe we know how to 
treat using asymptotic fields, will serve as a stepping stone to the nonabelian case. 
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